Limit shapes of Gibbs distributions on the set of integer partitions: 

The expansive case 



Michael M. Erlihson * 

Department of Mathematics, Technion-Israel Institute of Technology, 

Haifa, 32000, Israel. 

Boris L. Granovsky t 
Department of Mathematics, Technion-Israel Institute of Technology, 

Haifa, 32000, Israel. 

February 2, 2008 



American Mathematical Society 2000 subject classifications. 

P. '. Primary-60J27; secondary-60K35, 82C22, 82C26. 

> ■ 

•<-h ■ Keywords and phrases: Gibbs distributions on the set of integer partitions, Limit shapes, Random 
combinatorial structures and coagulation-fragmentation processes, Local and integral central limit 



theorems. 



* E-mail: maerlich@tx.technion.ac.il 

* E-mail: marl8aa@techunix.technion.ac.il 



1 



Limit shapes 



2 



Abstract 

We find limit shapes for a family of multiplicative measures on the set of partitions, induced by 
exponential generating functions with expansive parameters, ~ Ck p ~ l , k — > oo, p > 0, where 
C is a positive constant. The measures considered arc associated with the generalized Maxwcll- 
Boltzmann models in statistical mechanics, reversible coagulation-fragmentation processes and 
combinatorial structures, known as assemblies. We prove a central limit theorem for fluctuations 
of a properly scaled partition chosen randomly according to the above measure, from its limit 
shape. We demonstrate that when the component size passes beyond the threshold value, the 
independence of numbers of components transforms into their conditional independence (given 
their masses). Among other things, the paper also discusses, in a general setting, the interplay 
between limit shape, threshold and gelation. 

1 Introduction and Summary 

Given a sequence of probability measures {hn, N > 1} on the sets of unordered partitions of 
integers N > 1, the limit shape, provided it exists, defines the limiting structure, as N — > oo of 
properly scaled partitions chosen randomly according to the above sequence of measures. The study 
of the asymptotic structure of random partitions is stimulated by applications to combinatorics, 
statistical mechanics, stochastic processes, etc. Our paper focuses on limit shapes for a class of mea- 
sures associated with the generalized Maxwell-Boltzmann models in statistical mechanics, reversible 
coagulation-fragmentation processes and combinatorial structures called assemblies. In the course of 
the asymptotic analysis of the above class of measures if at, we reveal some interesting phenomenon 
that, as we expect, will be also seen in models related to other multiplicative measures on the set of 
partitions. 

We describe now the context of the present paper. In Section 2 we deal with an arbitrary 
probability measure on the set of partitions. We study here the linkage between the following 
three important concepts in statistical mechanics and combinatorics: threshold, gelation and limit 
shape. As a by-product of this study, we establish the non-existence of limit shapes for some known 
models. Section 3 gives a definition of a multiplicative measure and outlines the four main fields of 
applications of these measures. Section 4 contains the statements of our main results that are related 
to multiplicative measures induced by exponential generating functions with expansive parameters, 

Ofc ~ Ck p ~ l , k — ► oo, p > 0, where C is a positive constant. Namely, Theorems 4.1 and 4.3 

i 

determine the asymptotics of component counts of sizes that are of order NP+ 1 (=the threshold 
value), while Theorem 4.6 accomplishes the same for components of small sizes (=o (^Np+ t ^). In 
particular, in Theorem 4.3 we obtain limit shapes of Young diagrams under the measures considered. 
As a corollary of the above three theorems, we reveal that when the component size passes beyond 
the above threshold value, the asymptotic independence of the component counts transforms into 
their conditional independence, specified in Theorem 4.1. Section 5 provides proofs, that are based 
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on a far reaching generalization of Khitchine's probabilistic method. In the last Section 6 we discuss 
the limit shapes for our models versus the ones obtained by Vershik for the generalized models 
of Bose-Einstein and Fermi-Dirac. The two latter models correspond to a class of multiplicative 
measures induced by Euler type generating functions. 

2 The interplay between Limit shape, Threshold and Gelation 

We shall work with the set Qn = {v} °f an unordered partitions rj = (m, . . . , tin) '■ Ylk=i ^ n k = N, 
of an integer N. Here n/% = rik(rj) is the number of summands (=components) equal to k in a partition 
rj 6 S7jv- Each rj G can be depicted by its Young diagram (see e.g. pQ). The boundary of a Young 
diagram (shortly, Young diagram) of rj £ £In is a non-increasing step function v = v(»;rj) which is 
given by 

N 

(2.1) v(u) = u(u; rj) = } y n k (rj), u>0, rj £ Q N , 

k=u 

where we set rio(rj) = 0, rj £ £1^. By the above definition (|2.ip . for any integer u, the decrement 
v{u — 0) — v(u + 0) equals the number (possibly 0) of components of size u, z^(0) gives the total 
number of components, whereas the largest u with v(u) > equals the size of the largest component. 
Obviously, v(u)du = N ', r\ G Oat. Let now r = rjq > and for a given partition rj £ CIn define 
the scaled Young diagram V = V(»;rj) with the scaling (scaling factor) : 

(2.2) v(u) = r -jt-v(r N u), u>0, rj£tt N . 

When a partition rj £ {In is chosen randomly according to a given probability measure on O^r, it 
is natural to ask if there exists a scaling = o(N), such that the random curve v(u) converges, as 
N — > do, to some nonrandom curve. Such a curve, if it exists, is called the limit shape (of a random 
Young diagram). To give a formal definition, we will need some more notations. Let C = {/(•)} be 
the space of nonnegative non-increasing functions / on [0, oo) with J °° l{u)du = 1. Clearly, v £ C, 
for all rj £ Q^r. We supply C with the topology of uniform convergence on compact sets in [0, oo). 
For a given r^r denote by p TN the mapping rj — > v of Qjy onto C. Given a probability measure /ijv on 
VLn-, the mapping p rN induces the measure p TN PN on £, {p rN PN){l) '■= MaK/VatO) ' e ^' 

In this section, we refer to pn as an arbitrary probability measure on f2/v- The definition of a 
limit shape given below follows the one by Vershik in [44j . 

Definition 2.1 ^4 continuous curve I £ C is called the limit shape of a random Young diagram w.r.t. 
a sequence of measures {pn, N > 1} on {0,n, N > 1} (£/ie Zimii shape of pn) under the scaling 
■ = o(N), N — > oo i/ i/ie sequence of measures {p rN pN,N > 1} on £ weakly converges, as 
N — ► oo, to i/ie defta measure which is concentrated on the curve I. 
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Clearly, the weak convergence of the sequence of measures {p rN fj,N, N > 1} in Definition 2.1 is 
equivalent to convergence in probability ofu(u), u>0 w.r.t. the sequence {hn,N > 1}. Namely, 
the continuous function Z(-) 6 £ is the limit shape of the measure fijsr under a scaling rjy if for any 
< a < b < oo and e > 0, 

(2.3) lim ^n{v 6 : sup |?(u) — /(u)| < e} = 1. 

A^ 00 ue[o,6] 

(|2.3p expresses the law of large numbers for the random process v(u), u > 0. 

Note, that if a measure has a limit shape I G £ under a scaling r^r, then for any c > the 
function cl(cu), u > is a limit shape of the same measure /xjy under the scaling crjy. 
We will say that a measure has no limit shape if there is no scaling rjy = o(N), that provides (|2.3p 
for some I G C. 

A sketch of the history of limit shapes. The evolution of shapes of random ensembles of 
particles, as the number of particles goes to infinity, was studied for a long time in a variety of 
applied fields: statistical mechanics (the Wulf construction for the formation of crystals, see [39], 
[38]). stochastic processes on lattices (the Richardson model, see [10 J), biology (growth of colonies) 
etc etc. A special study was concentrated on limit shapes for random structures on the set of 
partitions, in view of applications to statistical mechanics, combinatorics, representation theory and 
additive number systems. In 1977 two independent teams of researchers, Vershik & Kerov [42] and 
Shepp & Logan [29], derived the limit shape of a Young diagram w.r.t. the Plancherel measure. 
Following this seminal result, Pittel [35] found the limit shape of Young tableaux w.r.t. a uniform 
measure. Since the number of Young tableaux corresponding to a given partition (Young diagram) 
is known to be equal to the degree of the irreducible representation associated with the partition, the 
above uniform measure, as well as the Plancherel measure, is related to the hook formula. It should 
be mentioned that the research in this direction revealed also a deep linkage to the random matrix 
theory, which is now a rapidly growing subject (see [32]). Parallel to this line of research, Vershik 
|44| developed a general theory of limit shapes for a class of measures he called multiplicative and 
which are discussed in Section 3. These measures encompass a wide scope of models from statistical 
mechanics and combinatorics, but do not include the measures associated with the hook formula. 
The results on limit shapes of multiplicative measures obtained by Vershik and Yakubovich [43 j - 
[45j . [50] . [48] during the last decade concern measures induced by Euler type generating functions. 
(We refer to some details of this research in the course of the present paper). Extending these results, 
Romik [36] derived limit shapes for multiplicative measures corresponding to some restricted integer 
partitions. Note that that the limit shape of the uniform measure on the set of partitions (which 
is a multiplicative measure) was firstly obtained via a heuristic argument, by Temper ley [H]. A 
comprehensive study of this case was done by Pittel in [M] . 

In contrast, the multiplicative measures /xjv considered in our paper are associated with exponential 
generating function. 
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This being said, it should be stressed that the results on limit shapes of multiplicative measures 
were stimulated by the remarkable papers of Erdos, Turan, Szalay as well as by other researchers, on 
statistics related to integer partitions (for more details see [16], [35], [32], [22]). In the course of research 
on limit shapes of multiplicative measures, links of the subject to various fields of mathematics were 
revealed. In particular, recently the application of probabilistic methods to the study of logical limit 
laws was implemented in [TU] and [40J . This link is based on fundamental theorems of Compton that 
were extended and deepen by Burris and Bell [Z],[6]. The place of limit shapes in this latter field is 
not yet understood. 

It turns out that the existence of a limit shape is closely related to two other phenomena known 
in statistical mechanics, which are gelation and threshold. We recall 

Definition 2.2 Let qN(v) be the size of the largest component in a partition r] = (m, . . . , n^) G fijv, 
i.e. qN{f]) = max{z : rij > 0}. 

(i) We say that a measure /Utv exhibits gelation if for some < a < 1, 

(2.4) lim sup fi N {r] G tt N : qN(v) > aN } > °- 

(ii) A sequence qjy = o(N), N — > oo is called a threshold for the size of the largest component, under 
a measure fJ-N ( shortly, threshold of a measure ^n), if 

(2.5) lim fi N {r] G £In ■ lN{rf) < Qn} = 1, 

N-*oo 

while for any sequence vn, s.t. vn = o((Jn), N — > oo , 

(2.6) lim sup [In {r/ G Qn : qN{v) < v n} < 1- 

Remark 2.3 (i) In physics and chemistry, gelation is viewed as a formation of a gel which is a two- 
phase system consisting of a solid and a liquid in more solid form than a solution. In combinatorics, 
equivalent names for gelation are connectedness of components and appearance of a giant component. 
We notice that in \2J$ , in contrast to our definition l[2.4\ ), the definition of gelation ( =formation a 
giant component) requires that 

(2.7) lim fi N {7] £ Q N : q N {rj) > aN} = 1, 

for some a > 0. In \21$ , some sufficient conditions for the absence of gelation (in the above sense) 
are given for combinatorial models encompassed into Kolchin's generalized allocation scheme (see 

mi)- 

(ii) It is common (see e.g. \11$ ) to seek a threshold q~N, if it exists, in the form q^ = N@ , where 
the exponent (3 := inf{/? : /iA?(r/ G O^r : Qn(i]) ^ = 1}. 



Limit shapes 



6 



We write (L), (G), (T) to abbreviate the statements: "There exists a limit shape/ there exists gela- 
tion/ there exists threshold", respectively and write (§) to denote the negations of the above state- 
ments. 

Proposition 2.4 Let be a measure on Oat- Then 

(2.8) (G) (T), and (L) (G). 

Moreover, if fi^ has a limit shape under a scaling r^, then has a threshold q~N > O(rjv), N — ► oo. 
Proof: 

By the definition of gelation, 

(2.9) (G) O lim ^ {r? G ^tv : gjvfa) > eiV} = 0, 
for any e > 0. Therefore, 

(2.10) (G) lim fi N {77 € fi* : q N { V ) < e N N} = 1, 

for some sequence eAr — ► 0, N — > oo. This proves the existence of a threshold when there is no 
gelation. For the proof of the second implication in (|2.8p . our strategy will be to show that (L) 
implies that for large N the major part of the total mass N is partitioned among component sizes 
of order 0{r^). For this purpose, we employ the following argument. For given d > 1,U\ > 0, let 
uq = 0, Ui = iu\, i = 1, . . . , d be d + 1 equidistant nodes, and assume that I G C is a limit shape of 
fiN- Then 

u d r N -l d-1 u i+1 r N -l d-1 I u i+1 r N -l \ 

( 2 - 11 )]V ^ kn k(v) = Jj^2 X fen *(»7) ^ I ^ X "AC 7 ?)]' ??G^Ar. 



Next, (g3J> gives for all 1 < i < rf - 1 

{iti+irjv-l "j 
7] £tt N : Uij^- ^ n k(v) > u i{j-i u i) ~ K u i+l) ~ e i) \ = 1) 
k= Ui r N J 

for arbitrary ej > 0, 1 < i < d — 1. Consequently, substituting Uj = wl, 1 < i < d, gives 



(2.13) lim /j N < 

TV— +00 



r] G £1tv : — ^ ^(77) > ^ (Z(uj) - l(u i+ i))uj - ui y^ iei ( = 1 - 



Now we write 

d-1 d-1 

^2 ( l ( Ui ^ ~ l ( u i+l)) u i = ^2( U i~ Ui-l)l(Ui) - l(u d )u d -l = 
i=l i=l 
d-1 



(2.14) = (ui+l ~ Ui)l(ui) - l(ud)u d -i > / l(t)dt - l{u d )u d _ 

i=l 



1- 
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where the last equation is due to the fact that the points Ui, i = 0, . . . ,d are equidistant. Since 
J °° l(t)dt = 1 and the function I is non-increasing and continuous on [0,oo), we have that for a 
sufficiently large d and a sufficiently small u\ > 0, 

(2.15) / l(t)dt > 1 - -, l(u d )u d -i < -, 

for any e > 0. Now it is left to couple u\ with the above ei, 1 < i < d — 1 and e, by setting, say, 
6j = i -3 and < iti < ^(Ya^i ^ 2 ) _1 > to conclude from (|2.13p that for any e > 



(2.16) lim /iTv < 



u d r N -l 1 
eQN: N ^ fe (r?) > 1 - e ^ = 1. 



k=uir?f 

Since r^r = o(N), N — > oo, this immediately implies G and, therefore, T. We also derive from (|2.16p . 

(2.17) lim jUtv {r? G ■ Qn^) < u x r N } = 0, 

N^oo 

which completes the proof. ■ 

The expression (|2.16p exposes the following meaning of the scaling factor r^r, that is not imme- 
diately seen from the law of large numbers (12. 3p . When N is large, almost all mass N is partitioned 
into component sizes in the range [uirjv, Ud r N], for a sufficiently small u\ > and a sufficiently large 
Ud < oo. Consequently, the boundary of the major part of the random Young diagram scaled by rjy 
acquires, as N — > oo, the shape of a nonrandom curve I. 

Proposition 2.4 is applied in Corollary 3.1, in the next section to prove the non-existence of limit 
shapes for certain classes of measures on the set fi_/v. 

3 Gibbs distributions and multiplicative measures 

3.1 Mathematical setting 

The subject of the present paper will be the following class of measures fipj on the set of integer 
partitions of N : 

a ni a™ 2 a nN 

(3.18) ns{r]) = (cjv)" 1 -. V" JV I ' V = (ni,...,n N )en N . 

nv.n2i ■ ■ ■ tin'. 

Here a = {a^} is a positive function on the set of integers, called a parameter function, and cm = 
c/v(ai, . . . , a/v) is the partition function of the measure ' 

(3.19) co = l, c N =Y: W » N >1. 

j-f n 1 \n 2 \...n N \ 

Following Pitman (see [S],[33J), we call the measures of the form (|3.18|) . (|3.19|) Gibbs distributions. 
Gibbs distributions are incorporated into the following general construction formulated by Vershik 
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who was motivated by applications to statistical mechanics. Let = Uat>i be the set of 
all integer partitions r\ = (ni, . . . ,nk, ■ ■ ■) and let s = {sk}f be a sequence of positive functions on 
the set of integers. We associate with a partition rj £ Qn the function Fj^(rj) = Y\k=i s k{ n k) and 
define the measure /xjy on by fijsfiv) = ( c n) Fn(v)i where cat = ^2 v£ q n Fn(t]) is the partition 
function of hn. Next, the family of probability measures (i( x > on f2, depending on a parameter 
x > is constructed in such a way that for each N > 1 the conditional probability of /j,( x > given 
Otv, is the aforementioned measure /xjv> f° r all x > from the domain of definition of the function 
!F(x) = X]jv>i cnx n . Explicitly, fi^ x > is given by 

^\r l ) = x N W(F(x))- 1 F N(Tl) (r l ), r, e 0, < x : ^(x) < oo, 

where N(rj) = Y2k>i ^ nfc denotes the number N which is partitioned by an jj G and T is the 
partition function of fi^ x \ or equivalently the generating function for the sequence {ctv}i°- It is not 
hard to verify that the measure possesses the required conditioning property which formally 
reads as follows: 

(3.20) m ( v ) := (^\n N )( V ) = ryGOiv, 

for all x in the domain of definition of the function T . Moreover, it follows from the above definitions 
that the generating function T is expressed as the Cauchy product of the generating functions 
^k( x ) = J2r>o s fe( r ) a - r > k > 1 for the sequences {sfc(r)} r >Q, k > 1 : 



(3.21) F(x) = l[F k (x k ). 

k>l 

Consequently, 

(3.22) ^ x Hv)=x N{v) \\ k=1 J) J , r/Gfi, < x : T{x) < oo 
and 

N 

(3.23) ^JV(^) = (cat) -1 Y\ s k(nk), V G 

k=l 

In view of the representation (I3.22h . Vershik calls the family (with respect to x) of measures y^ x > 
multiplicative. Following [35], we will preserve the same name for measures yUjv induced by multi- 
plicative families of measures yS x ' , via (I3.20p . The representation (13.221) tells us the important fact 
that the random counts nk, k > 1 are independent with respect to the probability product measures 

From the above formulae the following relation between the two families of measures y^ x ' and 
fj,N holds: 
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(3.24) = x N {^(x))- 1 c N /x N , 0<x: F{x) < oo. 

N>1 

This says that measures can be viewed as a poissonization of measures which is a standard 
way to deduce the equivalence, as N —* oo, of canonical and microcanonical ensembles in statistical 
mechanics (see [9], [Uj and references therein). 

In the above setting, the following three particular forms of the functions s k , k > 1 and corre- 
sponding to them generating functions T are of great interest: 

Case 1. 

a r 

Sk{r) = -7, r > 0, a k > 0, k > 1, 

(3.25) J^(x) = exp(^a fe x fc ), < 2 : ^ a fc x fc < 00. 

fe>l fc>i 

By (Ml, 

(3.26) ^ (77) = J] ( (afcXfc , )Wfc exp (-a fc x fc )) , 77 = (m, . . . , n k , . . .) € Q, 

k>i Uk ' 

i.e. is the probability product measure induced by the sequence of Poisson (akX k , k > 1) random 
variables. It follows from (|3.23p that in this case the measures hn, N > 1 are the Gibbs distributions 
defined by (|3.18p . (|3.19p . The associated generating function T is called exponential. 
Case 2. 



/rrifc + r — 1\ 

s fc (r) = ( ^ \ , r > 0, m k > 1, k > 1, 



(3.27) ^(x) = II fl _ Ifcw . 0<x:^m fe x fc <oo. 

k>l ^ ' k>l 

Consequently, (ET221) . (T3T23T) give 

^GOfo) = X N M Y[(l-x k ) mk , i)£0, 
fe>l 

(3.28) ^(r?) = ( CJV )- 1 TT ( mh + Uk " 1 Y £ Ojv, 

ii V n k J 



which says that p^ x > is the probability product measure induced by the negative binomial random 
variables NB(mfc,x fc ), k > 1, with a free parameter < x < 1. The function J 7 is called the Euler 
type generating function, because in the case m k = 1, fc > 1 it conforms to the standard Euler 
generating function for integer partitions. Note that in the latter case fijsr is the uniform measure on 
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Case 3. 



s fc (r) = II, <r < m k , m k >\, k > 1, 

(3.29) F{x) = \\{l + x k ) m \ < x : ^m k x k < oo. 

fc>i fc>i 

In view of (13.29[) . the following notation for the sets of restricted partitions is needed: 

fi(mi, . . . ,mk, ■ ■ •) = {f? = (^l, • • • ,«fc, • • •) G : < n fc < m fc , k > 1}, 

(3.30) Ojv(mi, . . . , mjv) = {?7 = (^i, . . . , rijv) € Ojv : < nfc < mfe, 1 < A; < N}. 

We have from (13321) . (13331) 

Uk>i( l + X ) k 



N ( \ 

(3.31) m ( V ) = (cat)- 1 !] m " , r/G^ 

fc=l Vnfc ^ 



(mi, . . .,m N ). 



Thus, is the product measure induced by the binomial random variables Bm(m k ,x k ), k > 1, 
with a free parameter < x < 1. In particular, if = 1, fc > 1, then 0(1, . . . , 1, . . .), 0^(1, . . . , 1) 
are the sets of integer partitions with distinct components, i.e. each n k is either or 1, and /i^v is 
the uniform measure on Sljv(l, ... ,1). 

We see that in all the aforementioned three cases the sequence {/ijv, N > 1} of multiplicative 
measures is given by a single parameter function, which is a = {a k } in the first case and m = {ink} 
in the two other cases. It is known that for a great scope (but not all!) of applied models, the 
asymptotics of a parameter function can be described by a power law with the exponent (p— 1) G H, 
by which we mean that the parameter function behaves asymptotically as Ck p ~ l , C > 0, k — > oo. 
From the point of view of the asymptotics of the corresponding models, the following three ranges 
of p should be distinguished (see [3]): the logarithmic case, p = 0, the expansive case, p > 0, and 
the convergent case, p < 0. It was found (see [3] and 0] respectively) that in the logarithmic and 
the convergent cases the measures /xjy exhibit gelation. In view of this, Proposition 2.4 leads to the 
following 

Corollary 3.1 In the logarithmic and the convergent cases, multiplicative measures have not limit 
shapes. 



Therefore, multiplicative measures may have limit shapes in the expansive case only. 
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3.2 Applications 

We outline the four main fields of applications of multiplicative measures. 

• Coagulation- fragmentation processes (CFP's). Given an integer N, a CFP is a continuous- 
time Markov chain on the set = {r]} of all partitions of N. Here N codes the total population 
of indistinguishable particles partitioned into nj groups of size j, j = 1, . . . , N. The possible 
infinitesimal (in time) transitions are coagulation of two groups (=clusters) of sizes i and j into 
one group of size i + j and fragmentation of one group of size i + j into two groups of sizes i 
and j. If the ratio of these rates is of the form 

(3.32) ^±i, 2<i + j<N, 

CLidj 

where a = {cij} is a positive function on the set of integers, called a parameter function of a CFP, 
then the process conforms to the classic mean-field reversible model of clustering formulated 
in 1970s by Kelly and Whittle. (For more details and history see [25], 09], p], [H], P2],[II])- 
It is known (see [25j,[llJ) that, given a parameter function a, the invariant measure /in of the 
corresponding CFP is exactly the Gibbs distribution (|3.18|) . In particular, under the parameter 
function at- = k > 1 among the possible rates of coagulation and fragmentation are the 

ones given by the kernels K(i,j) = {ij) p ~ l and F(i,j) = (i + j) p ~ 1 respectively. The cases 
p = 1, constant rates, and p = 2 are common in chemical applications (see [11], |22j). 

• Decomposable random combinatorial structures. A decomposable structure of size N 
is a union of indecomposable components, so that the counts m, . . . , rtjv of components of sizes 
1, . . . , N, respectively, form an integer partition of N. Given a sequence of integers m = {m^}, 
it is assumed that each component of size k belongs to one of the m/% types. Three classes of 
decomposable structures: assemblies, multisets and selections, encompass the whole universe 
of classical combinatorial objects. Assemblies are structures composed of labelled elements, 
so that afc = , k > 1, multisets are formed from unlabelled elements and selections are 
multisets with distinct components. Supposing that a structure is chosen randomly from the 
finite set of a certain class of structures with size N, the random partition of an integer 
N, called component spectrum, is induced: 

N 

K (N) = {K (N)^ .^ K (N) ): J2 kKi N) = N, N > 1, 

k=l 

where the random variable if! represents the number of components of size k in the random 
structure. The distribution of is described by a probability measure on the set of 

integer partitions of N. A remarkable fact is (see[3]) that the distributions of K.( N ^ correspond- 
ing to assemblies, multisets and selections are just the multiplicative measures fijsi induced by 
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Poisson, negative binomial and binomial distributions respectively. Based on the combinato- 
rial context, it is also common to treat distributions hn of as the ones generated by 
the conditioning relation (see [3], [18]), a frame which is close (but not identical!) to Vershik's 
formalism. 

• Statistical mechanics. In this context, the multiplicative measures /x x ' and are referred to 
as macrocanonical and microcanonical ensembles (of particles) respectively. For more details 
on this topic, see [Hj, [IB]. The multiplicative measures induced by Poisson, negative binomial 
and binomial distributions, with constant parameter functions, provide a mathematical setting 
for the three classical models of ideal gas, called Maxwell-Boltzmann (MB), Bose-Einstein 
(BE) and Fermi-Dirac (FD) statistics, respectively. In [H] the multiplicative measures with 
parameter functions Ck p ~ l , C > 0, p > 1, k > 1, are called generalized classical statistics. 

• CFP's on set partitions (see [8].[33j). This comparatively new setting is based on Pitman's 
study of combinatorial models of random set partitions, which developed from the Ewens 
sampling formula and Kingman's coalescence processes. We assume here that in the preceding 
setup for CFP's, particles are labelled by 1,. . . ,N, so that the state space of a CFP is the 
set fiyv] = {vr^j} of all partitions 7qjv] of the set [N] = {1, . . . , N}. Denoting \Aj\ the size of 
a cluster (block) Aj ^ [N] of 7T[jv]; we assign to each Aj, a weight m\A_.\ which depicts the 
number of possible inner states of Aj, the states can be e.g., shapes (in the plane or in space), 
colors, energy levels, etc. This says that to a set partition 7Tr n i jfc into k given clusters A\, . . . , A k 
correspond n^=i m |A,| different states of the CFP considered, so that the total number of 
states formed by all partitions of the set [N] into k clusters is equal to 

k 

(3.33) Y[m\Aj\ '= B N,k, 

where Onyifc is the set of all partitions of the set [N] into k blocks and the number B^^ is 
known as a Bell polynomial in weights mi, . . . ,m7v-fc+i- Next, for a given k > 1 the uniform 
measure P[N],k on the set On\nfc = {^[N]^} 1S defined: 

rij=i m \Aj\ 

(3-34) P[jV],fc(7r[iV],A:) = ^ —j ^[N],k & $}[N},k- 

In a more general setting which encompasses a variety of models (see [8j , [33! ) , the weights mj in 
(|3.34p are allowed to be arbitrary nonnegative numbers. By the known combinatorial relation, 
the distribution p\N],k 011 ^[N],k induces the distribution p^,k of cluster sizes \Aj\, j = 1, . . . , k 
on the set fijvfe °f integer partitions of N into k positive summands: 

N 

(3.35) PNjk ( v ) = (B Nik )- 1 l[(^f)^ r — v={ni,---,n N )en N>k , 
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where the partition function -Bjv,fc defined as in (|3.33p can be rewritten in the following form: 

n -. 

ven N:k j=i y i>' 

Pitman calls the measures PN,k, P[N],k microcanonical Gibbs distributions on the sets fi^fe and 
fJrjvifc respectively. The linkage of pN,k to the Gibbs distributions p^ in (|3.18p is expressed via 
the conditioning relation: 

PN,k = ^N\^N,k, 

where /^at is the Gibbs distribution given by (|4.37[) with aj = , j > 1. 

However, the generic model associated with set partitions of [N] involves a wealth of spe- 
cific problems (see [8]) arising from treating pn\n,fc! k = 1,...,N as marginal distributions at 
any time t, of irreversible time continuous markov processes of pure fragmentation (or pure 
coagulation) on the state space 

Finally note that Gnedin & Pitman [17] and Kerov |26j studied exchangeable Gibbs distribu- 
tions related to the Ewens sampling formula (the case p = 0). 

4 Statement of main results 

Our forthcoming asymptotic analysis is devoted exclusively to Gibbs distributions in the expansive 
case, i.e. when the parameter function a of the measure /xjv hi (|3.18p is of the form: 

(4.37) a k ~Ck p -\ C > 0, p > 0, k -» oo. 

We already mentioned in the previous section that this case conforms to the generalized MB statistics 
and to reversible CFP's with certain transition rates. To explain the relevance of the case considered 
to random combinatorial structures, recall that in the combinatorial context, Gibbs distributions 
with the parameter function (|4.37|) describe assemblies with the number of indecomposable compo- 
nents of size k, equal rrik = a k k\ ~ Ck p k\, C > 0, k — > oo. Examples of such structures with 
p = 1, 2 were given in [12J. Generally speaking, we can think that in an initial assembly having, say, 
rhk indecomposable components on k labelled items, each one of the rhk components is distinguished 
by one of bk types. For example, in permutations, each of (k— 1)! cycles on k labelled items is colored 
by one of colors; in forests, each of k k ~ 2 trees on k labelled vertices belongs to one of b^ species, 
etc. In this setting, the resulting assembly is composed of nik = bkrhk components of size k. Thus, 
under b^ = k a , where a > in the first example and a > 3/2 in the second example, the resulting 
assembly satisfies (|4.37p . 

Throughout the rest of the paper we assume that all random variables considered are induced by the 
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Gibbs distributions [in with the parameter function (|4.37p . 



We note that Vershik's approach [44j,[45j, [46J to asymptotic analysis of the generalized BE and 
FD models is based on the expression (|3.20p . which allows to replace the measure /ijv by the product 
measure pS- x \ with a properly chosen x — > 1. In contrast to it, our strategy is the straightforward 
study of the asymptotics of /xjy, as N — > do, using stratifications of the integer N and the total 
number of components. However, both approaches are based on common ideas rooted in statistical 
mechanics and the saddle point method. 

We first introduce some more notations. Consider two sets of random variables Vj and Kj, j = 
0, . . . , q > 1 that determine for a given r\ £ Oj\r stratifications of the total number of components 
v = u(rf) = J2k=i n k( 7 l) an d of the total mass N = J2 k =i ^ n k(v)i respectively. Namely, for given 
integers q > 1 and = Mq < Mi < . . . < M q < N < M q+ \ = N + 1 and a given ij £ fijv, we set 

Ah-l N 

u = u (Mi,r]):= n k (rj), Vj = Vj(Mj, 77) := n k (r)), 

k=M k=Mj 
Mi-l N 

K = KoiM^ri) := kn k (ri) = N - K lt Kj = K^M^r,) := £ kn k ( V ), 

k=M k=Mj 

(4.38) j = l,-..,q 
and denote 

(4.39) D = (v ,...,v g ), K = (K 1 ,...,K q ). 

We will refer to Mj, j = l,...,q as stratification points. The random vectors ~u , K induce the 
decompositions 1?* and K* of the random variables v and N — Kq respectively, into sums of q + 1 
and q disjoint parts: 

~v* = (i/q,...,u*) := (fo,^i - vi,... - V q ,Uq) , 

(4.40) i?* = (K!,...,Kl):={K l -K 2 ,...,K q „ l -K q ,K q ). 

We also set JTq = ifp- Clearly, 

Afj+i-l Aff+l-l g 

(4.41) !/;= ^ n fe , = ^ kn k , j = 0,...,q, ^K* = N. 

k=Mj k=Mj j=0 

Formally, ~~v* and K* can be viewed as the 1-1 linear transformations B : TR q+1 — > H <?+1 and 
£>i : 1R5 — > H 9 of the random vectors ~~v and -fT respectively: 

(4.42) ~v* = B~u, K* = B\K. 
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It is known from p3] that in the expansive case (|4.37p . the measure /ujv has a threshold N? +1 , p > 

and that the major part of the total mass N, when N is large, is concentrated on components of 
i 

sizes of order Np +1 . Therefore, in view of the note after Proposition 2.4, we take throughout the 
rest of the paper 

(4.43) r N = h^N^TT, h = (CT(p+i))& 

as the scaling factor for the limit shape in question. Here the constant h is chosen in order to simplify 
the expression for the limit shape. 

To establish the desired limit theorems, it is required to define proper scalings of the two pairs 
of random vectors in (|4.39p and (|4.40p . We denote by • the corresponding scaled quantities. The 
explicit expressions for • will be chosen in accordance with the asymptotic problem considered. 
Firstly, in the forthcoming Theorems 4.1 and 4.3 and Corollary 4.5, we will study the asymptotics 
of the random quantities in (|4.38p . (|4.4ip in the vicinity of the rj\r. In this case, we consider the 
stratifications induced by the points Mj = [ujriy], j = 0, . . . ,q, M q+ \ = N + 1, where = u$ < 
u\ < . . . < u q do not depend on N, and take the scalings • in the form 

f>j = v7o(p-i) 



"1 








V 


% ' 




-El 










V 3 ' 


-Si 




^ ' 







(4.44) K* = VMp + T) J , 3=0,..., q 
and 

% = Vfo(p- 1) 

(4.45) Kj = V/o(p+l) JVJ ^ J \ j = 0,...,q. 
Here the constants 

(4.46) fj (s) = C x s e- x dx, j = 0,...,q, s>0, 

JUj 

whereas the quantities Sj,Ej,Vj depend on the parameter 5 = 5n = (^Af) _1 : 

A/i-l Mi-1 Mi-1 

S = Yl a k e ~ 5k ' E 0= Yj ka k e ~ Sk > V 0= Yl 

k=M k=Mo k=M 

N N N 

Sj = Y a * e ~ 5k > E i = E ka ^ 5k , Vj = Y k 2 a k e- 5k , 



■6k 



(4.47) 



j = 1, . . . q. 
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Finally, the starred quantities S*, E* in the RHS's of (|4.44p are defined respectively as the transfor- 
mations B and B\ of the corresponding non-starred quantities Sj, Ej. Also, we will use the following 
abbreviation for the conditional probabilities 



P W : 



•2,J 



when the particular values •i,j,»2,j are not important. 

The theorem below that deals with the asymptotic behaviour of the scaled quantities defined by 
(|4.44p . will be a source of our subsequent results on limit shapes. 

Theorem 4.1 Under N — > oo and a given q > 1, 

(i) The random variables 0*, j = 0, . . . ,q are conditionally independent given KJ,j = 0, . . . ,q, 
s.t Td=o K j = N - Moreover, 



(4.48) 



P (P\T*)=f[p(»*\K] 



(ii) The (2q + 1)— dimensional random vector {u*,K*) = (z>q , . . . , £>*, K^, . . . , K*) converges 
weakly to the (2q + 1)— dimensional gaussian random vector with zero mean and the covariance 
matrix Tn,q = 0, . . . , 2q} given by 



fm(p - l)lr 



, _ f m (p)h(p) 

=* r(p+2) ' 



(4.49) r mk (g) = { 



fm(p)^-m=k-q T(p+2) ' 

fl-Mlu fk(p)fm- q ( P +l) 

J k\l> ) s -k=m-q F(p+2) ' 

f r«-i-i"\i f m - q {p+i)fk- q (p+i) 

K Jm—q\F~T ± J J -m=k T(p+2) 



if < m, k < q, 

if 0<m<q, q + l<k<2q, 

if < k < q, q + \<m<2q, 

if q + 1 < m, k < 2q, 



where li = j is the Kroneker symbol. 

(hi) Moreover, moments of the random vector (i)*,K*) converge to the corresponding moments 
of the gaussian random vector in (ii) . In particular, 



(4.50) 



#mfc(?) 



lim < 

N^oo 



Covty^vfy, if0<m,k<q, 
Cov(v* m , K*_ q ), if 0<m<q, q + l<k<2q, 
[Cov(K^_ q ,K*_ q ), if q+ 1 < k,m < 2q. 



(iv) Let the stratifications be induced by the equidistant points Uj, j = 0, . . . ,q and let k = 
m + s, s > 0, < m, k < 2q. Then the absolute value of the covariance, l^^fc^)!' m onotonically 
decreases in s. In particular, if q — > oo, s — > oo, while m is fixed, then 



(4.51) 



O (u£exp(-u fc )) , 
O (u p k + _\^v{-Uk- q , 



if < m < k < q, 
if q < k < 2q, s 



oo. 
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Remark 4.2 The first part of h4-51\ ) expresses the exponential decay of correlations between the 
scaled counts of components of different sizes, as the "distance" s between the sizes goes to infinity. 
Phenomena of such kind are widely known in equilibrium statistical mechanics (see e.a.]21^. fWjj and 
references therein). 

Let I £ C be a limit shape of a measure hn and let v{u) be denned as in (|2.2p . We call 

(4.52) A(u) := D(u) - l(u), u>0 

the random fluctuation of //jy from its limit shape at a point u. To state the next theorem, let 

(4.53) b r (u) = CT(r + l,u), u>0, r > -1, 

where F(r + 1, it) = f e~ x x r dx, u>0, r > — 1 is the incomplete Gamma function. 

Theorem 4.3 (Limit shape of // v £fte cental limit theorem for scaled fluctuations from the limit 
shape). Under N —* oo and a given q > 1, 

(i) XTte random vector (Oi, . . . ,u q ) defined by ft4-45\ ) weakly converges to the q— dimensional 
gaussian random vector with zero mean and the covariance matrix {e m k}\ given by 

(4.54) e mfe = ^(u.) - m? fc = x 

r(p + 2) 

where s = max(k,m). 

(ii) Setting l(u) = Y(p+i) > u — ^ * n i4-52\ ), the relation between the scaled quantities i>i,...,v q 



and the random fluctuations A(itj) is given by 

(4.55) %~ (CT^ + l)^ 1 ) 1 ^^), j = l,...,g, iV ^ oo, 

/or all 7] £ {In. 

(iii) TTte measure has the limit shape 

r(p + i) 



?/■» 



der i/ie scaling rjy given by H4-439 - 



Remark 4.4 (a) Recall that the measure [in considered has the threshold Np+t . By (ii) of the above 
theorem and fl^.^3| ), i/te threshold turns out to be of the same order as the scaling rjy for the limit 
shape. In view of this, we believe that the following stronger form of the second part of \2. 8\) in 
Proposition 2.4 is valid: For a wide class of multiplicative measures, the existence of a threshold qjy 
implies the existence of a limit shape under a scaling rjy = 0(qjsr), N — > oo. 

(b) In \49j , an analog of our cental limit theorem for fluctuations was established for the uniform 
measure on the set of unordered partitions with distinct summands (=classical FD statistics). Recall 
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that this multiplicative measure is associated with the generating function Y\k>i(^ + x )~ ■ The 
structure of the covariance matrix in \43j is similar to the one in ft4-°~4\ )- In particular, the exponential 
decay of covariances is also seen there. 

Corollary 4.5 (The functional central limit theorem for the scaled number of components). Let in 
the above stratification scheme, q = X, u\ = u and denote v{u) = v\. Then, for all u > 0, the scaled 
random variable D(u) weakly converges to iV(0,en), where en = en(u), u > is as in \4-$4\ )- 

Recalling that z^(0) is equal to the total number of components in a random partition, the particular 

case u = of the above corollary recovers the central limit theorem for i/(0) established in |12j . 

We now turn to the asymptotic behavior of the counts of components of sizes o(qtv) (=small sizes 

i 

in comparison to the threshold), s.t. o((Jn) — > oo, N — ► oo, where q~N = Np +1 . For this problem we 
make use of the stratification points Mj = Mj(N), s.t. 

Mj = o [Np+i ) oo, N oo, Jim - ^-:= Pj <l, j = l,...,q, 
(4.56) M = 0, M q+1 =N + 1, 

while the scaled quantities £>*, K* are taken in the form slightly different from the one in ([4.440 : 



(4.57) 




Theorem 4.6 (Asymptotic independence and the central limit theorem for the counts of components 
of small sizes). Let the stratification points satisfy \4-50( ) and C>* = (z>q, . . . , v*), K* = (K*, ■ ■ ■ , Kg) 
be given by |^.57[ ). Then, as N — > oo, 

(i) The coordinates of each one of the two random vectors v* and K* are independent random 
variables. 

(ii) The random vector ~x = (Oq, Kq,0^, . . . K* 1 ,0*) weakly converges to the (2q + 1) 
-dimensional gaussian vector with zero mean and the covariance matrix having a diagonal block 
structure, with q blocks 



(4.58) BJ 1 :=[ 3 , j = 0,...,q-l, B~ l = 1 - a% 
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where 



3 ^ r J' 

(4.59) al = 1 - lim Var{P*) - M/;+l) - v ~ 



9 iV->oc v qJ T(p)T(p + 2) p 2 +p 

Remark 4.7 It was proven in \1J$ that in the model considered, the counts of components of fixed 

sizes (= the random variables nfe 1 ,...,rtfe,, < k\ < . . . < ki < oo, I > 1) are asymptotically 

independent. Combining this result with our Theorems 4-1 an d 4-6 says that when the component size 

i 

passes beyond the threshold value Np +1 , the asymptotic independence of component counts transforms 
into their conditional independence (given masses). In this respect, the threshold value can be also 
viewed as the critical value for the independence of component counts in the model considered. 



5 Proofs 

5.1 Khintchine-type representation formula 

As in [12], [13]) 03] an d [IS]) our tool for the asymptotic problems considered will be the probabilistic 
method by A.Khintchine, that was introduced in 1950's in his book [27]. The idea of the method is 
to construct the representation of the quantity of interest via the probability function of a sum of 
independent integer valued random variables depending on a free parameter, and then to implement a 
local limit theorem. In the course of time, this method, sometimes without mentioning Khintchine's 
name, was applied to investigation of a large scope of asymptotic problems arising in statistical 
mechanics and in enumeration combinatorics (see e.g. [28J , [2\ , [12] , [13] , [H] , [15] , [19] , [31] , [30] and 
references therein). The first three subsections of the present section contain preparatory asymptotic 
analysis towards the proof of our main results. 

Clearly, the distributions of the random vectors ~u and K are completely determined by the measure 
[AN . Explicitly, for any given vectors L = (Lq, . . . , L q ) and iV = {N\, . . . , N q ), 

P(7? = L) =^>(V = ~L,K = N y 

N 

(5.60) = J^P (p = = iV*) := J^i2(Z*,]v*), 

where, in accordance with (|4.40p . (|4.42p . we denoted 

^ = (Lo, ■ ■ ■ , L^) = B L , 

(5.61) iV 5 = (Nl,...,N q k ) = B 1 N. 

(|4.40p implies that for any vectors L* and iV*, 

(5.62) R0, N*) = P (V = ~£, K = N 
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Consequently, it follows from (|3.18p . (|3.19p that 

/ N n k \ 

(5.63) R{L*,N*) = c] f 1 E 

r]en N \k=i k ' J 

where 

(5.64) A = Ait^N 1 ) = {r]£Q N : z?(t?) = £*, ^(r?) = N*}. 

Our first goal will be to derive the Khintchine type representation for the probability R(L*,N*). 

(i) 

Setting ao = 0, we construct the array of discrete random variables j3} defined by 

(5.65) P($ ] = k) = ^f^, k = M j ,...,M j+1 -l, j = 0,...,q, 1>1, 

where the stratification points Mq, . . . , M q+ \ are as in (|4.38p . 5q > 0, . . . , 5 q > are free parameters 
and 

M j+1 -1 

(5.66) S* = S](Sj)= a ^~ &3 ^ J = 0, ■■-,<!, 

k=Mj 

in accordance with (|4.47|) . We assume that for a given < j < q, the random variables ft , I > 1 

(i) 

are i.i.d. and for a given / > 1, the random variables ft , < j < q are independent. 

Distributions of the type (|5.65p are widely used in the asymptotic analysis related to combinatorial 
structures, [28] , [3] , [30] . These distributions firstly appeared in the papers of Goncharov (1944) (for 
references see |28j ) and in the monograph [27] by Khintchine (1950). 

Lemma 5.1 (Khintchine 's type representation for the probabilities R(L*,N*)). Denote 



(5.67) Y®=Y,Pi' j = <>, ■■■,<!> 

i=i 

where f3\ j \ I > 1 are i.i.d. random variables given by \5. 65\) , [5. 66\) . Then 



-» — -> i [S* 

(5.68) R{L*,N*) = c„ 1 ]J± 7 ^—ex V (6 j N*)F(YW=N*), 



j=0 H 



where 



(5.69) N£ := N -^N*. 

3=1 
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Proof: It follows from (|5.65j) - ()5.67j) and the multinomial enumeration formula for the number of 
ways to distribute L* numbered balls over L*- numbered urns, that 



(5.70) 



3=0 



n p(f« = jv?) = n t^-p hw) I e (n^T 1 



iV 



i=o 5' 



ijEOn \k=l 



n k \ 



By ^MD , ([EMD this gives (^68|) . ■ 

Note that in contrast to the standard one parameter scheme of Khintchine's method, the represen- 
tation (|5.68|) is based on (q + 1) free parameters 5i, . . . ,5 q . 

Remark 5.2 One can see that \5. 68\) is a generalization of the representation formula for the total 
number of components which was obtained in 



5.2 The choice of free parameters 

Firstly, we set in (|5.68|) do = 6i = . . . = 6 q := 6 and apply the stratifications 

S ={S Q ,...,S q ), E = (E Q ,...,E q ), V = (V , . . . , V q ), 

as defined in (|Q7|> . of the quantities S := J2k=i a k e~ &k , E := J2 k= i ka k e- 5k and V := J2k=i k 2 a k e~ 5k 
respectively. Next, we take in (|5.60p the vectors L = [Lq, . . . , L q ) and N = (N\, . . . , iV 9 ) in the 
form 



(5.71) 



Lj = Sj + x jV r S , j = 0,...,q, 
Nj = Ej + Xj+gy/%, j = l,...,q, 



(5.72) 



Xj = < 



where Xj, j = 0, . . . , 2q are arbitrary reals. We also set Nq = N — Ni and adopt the following 
notation: 

xj, ifje{o,g}, 

Xj-xj+i, if 1 < i < g — 1, 

-x q+1 , if j = g+l, 

Xj-i — Xj, if q + 2 < j < 2q, 

x 2q , if j = 2q + 1. 

2q+l 

It follows from (|5TT2"j) that V] x° = 0. By ([BUT]) and (f5TTTj) we then have 

j=q+l 

(5.73) 
and 
(5.74) 



N* = E) + x° j+q+1 y/V , j = l,...,q. 



L* = S] + x°VS , j = 0,... 
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Note that (15391 . (ETTgp together with (Q71) imply that 



(5.75) 



N o= N ~Yl { E J + 



N-Et+xZ^y/Vo. 



We start by investigating the asymptotic behavior of the probability R(L*, N ) given by (|5.68p . when 
iV — > oo and 5 is fixed. We adopt from [13] the following representation (|5.76p of the partition function 
cat defined in (|3.19|) . Consider independent random variables s.t. k~ 1 ^ ~ Po(af 1 e~' 7k ), k > 1, 



v 



where <r > is a free parameter, and let Z^r = Then 



k=l 



(5.76) 



cjv = P(Zjv = iV)exp(iV<7 + S(<7)), N >1, 



N 



where S = S(a) = a^e ak = ^^5|(<r). We set <r = 5 . The forthcoming asymptotic analysis 
k=l j=o 

relies on the following fact that features the expansive case (|4.37p considered. Let 5 = 5n — > 0, as 
N —* oo and let the stratification points be as in Theorem 4.1, i.e. Mj = Ujr^, j = 0, . . . ,q, = 
Uo < Mi < . . . < u q . Then 



(5.77) 



lim S* = oo, S* = O(S ), N -» oo, j = 0,...,g. 



This is proved with the help of the integral test, as it is detailed in (|5.92p below. We will see from 
the proof of Theorem 4.6 that (|5.77p fails for the asymptotics of component counts of small sizes. 



Lemma 5.3 Suppose (5. 77 ) holds. Then for Lq, . . . , Lit given by {5.73 ) and for any 



N$,...,N*eN:J2N* = N, 



j=0 



(5.78) R(L*,N* 

where 
(5.79) 



(27T) 



- \j=o \ S j 



exp | -lj2( x °j) 2 ^ 



W 



nu F ( YU) = N t 



P(Z N = N) 



Proof: We have from (f5U8l) and (f5776l) 



rr 

i=o 



iV -> oo, 



exp (<5JV) exp ( - (5N + 5)) 



(5.80) 



JL (Sf) 
^exp(-S)n^ 



3=0 V ^ 
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Next, we apply Stirling's asymptotic formula to estimate (L*)l under (|5.73p and the condition (|5.7T|) : 



(Lp ~ J2*L) (L*)^exp(-L*) ~ J2nS* (S*p 1 + 



exp(-L*) 



27ns? {Sf) L ** exp (/.* log (l + | <<xp( - L) ] 



(5.81) 



- J2ttS] (S*p e~ b * exp 



N -> oo, j = 0, 



Substituting in (|5.80p the preceding asymptotic expansions yields ([5.78P 



We will follow the principle of Khintchine's method that the free parameter 5 should be coupled 
with iV so that to make the probabilities in the RHS of (|5.79|) large. Namely, (see also [13] , [14] ) we 
choose 5 = 5n as the solution of the equation: 

N 



(5.82) 



5k 



N. 



k=i 



We see from ([PT]) that (pT82j) implies E$ = E = N-E 1 . So, we obtain from ([5775]) and ([5772]) that 

(5.83) Nq = N — Ni = Eq + x° q+1 
which enables us to rewrite (|5.74p in a unified way: 

(5.84) AT* = E] + x° j+q+1 j = 0,...,q. 

It is easy to see that if a^, k > 1 are positive, the equation (|5.82p has a unique solution for any N > 1. 



Our next goal will be the establishment of local limit theorems needed for the proof of Theorem 



4.1. 



5.3 The local limit theorems for j — 0, . . . , q. 

Lemma 5.4 (The local limit theorems for j = 0, . . . , q). 

Let \Jt-.31 ) hold and let Mj = [ujrj^],j = 0, ... ,q, where = uo < u\ < . . . < u q , M q+ \ = N + 1 and 
[•] is the integer part of a number. Then for L*, N* as in ( 5. 73), |5. 74\ ) and for S given by ([57 

i / ^ 

(5.85) P(y (j) = N*) ~ 



1 exp f ] , N -> oo, j = 0, 



where 



lim 



are positive constants calculated in (5.107) below. 
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Proof: We denote by (pW the characteristic functions of the random variable to obtain 
(5.86) ¥{Y {j) = N*) = (2TT)- 1 ^\t)e- iN Pdt := (2tt)- 1 I^, j = 0,...,q. 



We will focus now on the asymptotics, as ./V — > oo, of the integrals 1^', j = 0, ... ,q. For any 
< \to t j\ < 7r the integral 1^) can be written in the form: 

(5.87) I^=l{ j) +4 J \ j = 0,...q, 

where I± = Ii\to,j) and 1% = I^i^oj) are integrals of the integrand of taken over the sets 
[— toj,toj] and [— it, —toj] \J[to,j,^], respectively. Using the technique of [13], [Hj, we will show that 
for an appropriate choice of < to,j = to(N, j) — * 0, N — > oo, j = 0, . . . , q, the main contribution to 
I®, as N — > oo, comes from Jp , i.e. from a specially constructed neighborhood of i = 0. Firstly, 
observe that 

(5.88) ^\t) = (y?(t)) Lj , j = 0,...,g, 

where ip^ is the characteristic function of the random variable (3j := /3^' defined by (|5.65p . (|5.66p . 

(i) 

To derive the asymptotics of the integral J) , under the above choice of to j, we will look for the 
approximation of ip^'(t), as t — > and N — > oo. For this purpose we need the asymptotic expressions, 
as iV ^ oo, for E/%, E/3] and E/3?. We have 

(5.89) E^ = -I, E/?J = ^, E/3? = -=(-, j = 0,...,g, 

^ ^ 



where ilj = H*{5, Mj, M j+1 ) = ^ k 3 a k e~ Sk , j = 0,...,q. It follows from the definitions of 

k=Mj 

Sj, Ej, and Hj that in the case considered the problem reduces to estimation of sums of the form 

Mj+i-l 

(5.90) a kk r e~ Sk , a k ~Ck p -\ k ^ oo, p>0, C > 0, r > 0, 



when N — > oo, <5 is given as the solution of (|5.82p and Mj, j = 0, . . . ,q + 1 are as in the statement of 
the Lemma. The asymptotic solution of (|5.82p was obtained in |14j : 

(5.91) 5 = 5 N ~ hN~^ = (rjv)" 1 , TV -» oo, 

where /i, r^r are as in (|4.43j) . Thus, Mjd — > Uj, N — > oo, j = 0,...,q. In view of (|4.37p . it is 
convenient to write a k = Ck p ~ 1 G(k), where the function G is s.t. lim G(k) = 1. Then, applying 

k—rOO 

the integral test (=Euler summation formula), we get from (|5.9ip 

V a k k r e~ Sk ~ C/ x p+r - 1 G{x)e- 5x dx = C5~ p ~ r x p+r - l G[^)e- x dx 

C5- p - r x p+r - l G(%\e~ x dx^5- p - r C x^'^e^dx 



<5, 

(5.92) := <5- p -7i(p + r--l), iV ^ oo, j = 0,...,(?, r > 0, p > 0, 
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where we set u q+ \ = oo. Note that by our notation (|4.53j) . 

(5.93) fj(s) = bsiuj) - b s (u j+ i), s > -1, j = 0,...,q. 
We also obtain from (15321) and (|5T89jh 

(5.94) E/?J = 0(r r ), iV^oo, r>0, j = 0,...,g 
and observe that 

(5.95) f J ( p + l)f j (p-l)-(f J (p)f > 0, / i (p + 2)/ i (p-l)-/ i (p + l)/ j (p)>0, i = 0,...,g. 

The first of these inequalities follows immediately from the Cauchy-Schwarz inequality, while the 
second can be derived by substituting fj(p — 1) from the first one and then applying again the 
Cauchy-Schwarz inequality. 

To arrive at the required asymptotic formula for ip^\ we first write for a fixed N, 

(5.96) <4 j) (t) = 1+iEPjt- ^E/3|t 2 + 0(E/?|t 3 ), t -► 0, j = 0,...,g, 
and then couple i with TV by setting 

(5.97) t 0J = 5(L*)-hog 2 5 = 0(5 1+P 2 log 2 6), j = 0,...,q, N -» oo, 
where the last equality follows from (|5.92|) and (|5.73p . Consequently, 

(5.98) (E/3J) t^. = 0(5^ log 2r 5) -»• 0, N -» oo, j = 0, . . . , 9 , r > 0. 
Now ([5T88]) and (15T961) together with (l5T73|) . (15T7H) give 

^(t)^* = + iE p jt _ I E /? 2 i 2 + 0(E^t 3 ))^e- iJ V* 

~ exp fti (LjEfy - NJ) - ^L*Bp]t 2 + i^t 2 (E^) 2 + L*0(E/3 3 i 3 )^ 

= exp ^it (EY (j) - TV/) - it 2 Vary (j) + L*0(E^ 3 )^ , 

(5.99) |*|<*o,j, iV^oo, j = 0,...,(?. 

In view of (|5,94p and (j5.95p , the choice (|5.97p of t j provides 

(5.100) lim t 2 Q , VarY^ = oo, j = 0,...,q 
and 

(5.101) lim ml ,E/? 3 = 0, j = 0, . . . ,g. 
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Remark 5.5 By fiTS® and /PS}) . 

(5.102) VarY® = L* (e/3 2 - (E/%) 2 ) = 0(<T 2 L*) , N -> oo, j = 0, . . . , g 
and 

e(f« - EF^) 3 = L*E(/5j - E/^) 3 = L* (e/? 3 - 3E/3? E^- + 2(E/5 j ) 3 ) = 0(<T 3 L*) , 

(5.103) iV^oo, j = 0,...,q. 

Hence, the following weaker ( =the third moment is not absolute ) form of Lyapunov 's sufficient con- 
dition for the convergence to the normal law holds for the sums j = 0, . . . , q: 

(5.104) Jim g^- = 0, j = 0, . . . , q. 

N ~*°° (Vary(?"))3 

This explains the existence of io,j that provides \5. 1 00}) and \5.101}) . The phenomenon described 
above is typical in applications of Khintchine's method (see \1S$ , \1J$ ,}1 9$ )- 

To continue the asymptotic expansion (|5.99p . we obtain from ()5.92j) . (|5.89p and (|5.95p 

(5.105) iV^oo, j = 0,. ..,<?, 
where 

f 2 (p) 

(5.106) = >0 , j = 0,...,q. 

The asymptotic relations below are the consequence of f|5. 105|) , (|5.73p and (|5.74p : 



] 3 J J o J _ 

X 4 , — Xj +q+1 - 



o /J(P) / /0(P-1) o / /0(P + 1) ._ f O O V 



7i(p-i) , 

(5.107) iV^oo, j = 0,...,g. 
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Now we are in a position to evaluate the integrals I± \ j = 0, . . . ,q. By virtue of (|5.99p - (|5.10ip and 
(I5T0711 . 



r(j) 



to, 



— *0, 



ip^\t)e- iN Pdt 

(EyO) - at*; 



exp 



y* ' 3 exp ^t(Ey^) - TV*) - Vvai-y^ 



2VaryO') 



exp 



0,J 



-- t\/Vary(j) 



(5.108) 



1 

V2¥ 



WaryO') 



exp 



exp 



exp 



(EY® - Nff 
2VaryO) 

(EY® - Nff 
2VaryO) 



•y/VaryO')t ,_ 
VVarYO)t j 



exp 



EyO') - at*' 

5 WaryO') 



-^(t-i^) 2 



exp 



-^(t-i/%) 2 



df 



N^oo, j = 0,...,q. 



Next we turn to the estimation, as ./V — > oo, of the integrals , j = 0, . . . , g. We start with 



r(j) I .. 



to. 



< 2 



(5.109) 



47T 



1,J 



Afj+x-1 



(27rit-5)fc 



k=M, 



f \cp( j \t)\dt = A7r ! 



"idt, j = 0,...,q, 



\<pV\2irt)\dt 



where we set tij = > 0, j = 0, . . . , q. Denote 

(5.110) g®(t) = \<p¥\2nt)\, t€ [t ld , 1/2], j = 0,. ..,</. 

It is easy to see from the definition of the random variables (3j that g^\l/2) < 1, j = 0, . . . , q. Since 
the (3j, j = 0, . . . , q have lattice distributions with span 1, this implies (see [10] . p. 131, [37], p. 286) that 
for a fixed iV, we have g (j \t) < 1, t G [tij, 1/2], = 0, ... ,q. Moreover, s> (j) (iij) 1, N -> oo, 
because iij — > 0, iV — > oo. We then conclude that for sufficiently large N and any < t < 1/2, 



(5.111) 

By the same argument, 



,0") 



(h d )>g®{t), j = 0,... 



0, 



(5.112) 



iV^oo, j = 0, 
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This and (|5.7ip enable us to write 

1 + ( 5 i) _1 2 a k e~ Sk (e 27Tiktl 'i - 1) 



// ; "i'i,: V 



/ Afi+i-l 

~ exp -2(1 + 0(55)) ^ a A e- 5fe sin 2 (ttHi,^ 
\ 

(5.113) iV^oo, j = 0,...,(?. 
Denote 

(5.114) D j = 2 a k e- Sk sin 2 (nkti d ) j = 0,...,q 

k=Mj 

and let Uj G (iij, Uj+i), j = 0, . . . , g. Then we have 
Dj > 2 a k e~ 5k sin 2 (ir kt ld ) > 2 min (a fe e _<5fc ) ^ sin 2 (vr^ij 

K=UjN P+ 1 v ' K=UjNP +i 



(5.115) 0(<T P+1 ) ^ sin 2 (7rHij), iV -» oo, j = 0,...,g. 

i 

k=ujNP+ 1 

Denote Cj = iij — Uj > 0, j = 0, . . . , q. To estimate the last sum in (|5.115p . we employ the following 
inequality from [15]: 

l+m—l ^ j 

(5.116) 2 ^ sin 2 (7rf/<0 > — min{l,(im) 2 }, |t| < -, m > 2, Z > 1. 

k=l 

i i 

We apply (|5.116p with I = ujNp+i and m = ejNp+ l + 1 to get from (|5.97p and (|5.9ip . 



1 ' 

2 V] sin 2 (7r/«tij) > -e^i^ min{l, (ejiV^itij) 2 } 
i 



(5.117) = O^-^log 4 ^ N^oo, j = 0,...,q. 
Now we deduce from (|5.115p . (|5.117|) . (I5.1U9|) and (15.1111) that 

(5.118) < 0(exp(-log 4 5)), iV -> oo. 
Comparing this with (|5.108p gives 

(5.119) I^ = (l{ j) ), iV^oc, j = 0,...q 
which together with (|5.107p . (|5.108p proves the lemma. ■ 
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5.4 Completion of the proofs of Theorems 4.1, 4.3 and 4.6 



To complete the asymptotic analysis of the probability R(L*, N*) in the vicinity of r^, it remains 
to derive the asymptotics of P(Zn = N), as N — > oo. It was found in [13] that 

1 



(5.120) 
where 
(5.121) 



¥(Z N = N) 



V2vrVarZ A r 



N oo, 



jv 



VarZiv = ^ k 2 a k e- 8k ~ r(p + 2)<T (p+2) , iV -» oo. 



fc=l 



Substituting in (15. 78ft the asymptotic expressions (|5.85p . (j5.105j) . (|5.120p and (j5.121j) we obtain 

<2 



R(L*, N ) 



y/VaFZ 



(2tt) 



9+ 



? n 



i 



3 -=o x /S|VarY(i) 



j(q(p+l)+p/2) V r (P + 2 ) | -Q 



(2tt)< 



exp 



cxp 



(5.122) 



iV ^ oo. 



Next we employ the scaled quantities z>*, K* constructed in (j4.44jl and denote 



(5.123) 



o ( ° °\ ° ( 

x — v^O' • • • ) x q)i y — \ x . 



9+2 > 



Note that K * = - Y? j= i Kt, by §EM). We get from (RTT22|> 



y^¥+2) f fo(p-l) \ + 19 



(2tt) 2 ^ 1 V So 



(5.124) exp 



" (^°) 2 
2^?/>-l) 



1 9 

4e 



i=o 



exp 



x 



iV ^ oo, 



^q+l — Sj=l X g+j+l- 



where x z q 

Each term of the sum in the second exponent depends on the two variables x° and only. 
Hence, under a fixed vector y° , the exponents in the RHS of (I5.124j) factorize into a product of q + 1 
terms, each one depending on only one of the Xj, j = 0, . . . , q. This proves the claim (i) of Theorem 
4.1. 

For the proof the claim (ii) of Theorem 4.1, we firstly need to deduce from (|5.124|) the central limit 
theorem for the (2q + 1)— dimensional random vector (u , K*). We see from the definition (|5.107p of 
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Kj, j = 0, . . . ,q that the expression in the product of the exponents in (|5.124p is a quadratic form of 
the coordinates of the vector (x° , y°). For a given q > 1, we denote by @*(q) the (2q + 1) x (2q + 1) 
matrix of this quadratic form. Let now Cj < dj, j = 0, . . . , 2q and 

(5.125) (c , d) = (c ,d ) x (a,di) X . . . X (C2q, d2q)- 

We also define the sets of discrete points 



(5.126) 



G 



U) 

N 



{z e {c^df) : S* + z^Jj^ G N}, if < j < q, 
{z E : Ei + zjj^ e N}, if q + 1 < j < 2q 



and Gn = G$ x . . . x G§ 9) . 



Lemma 5.6 ( The central limit theorem for the vector (£>*,K*)). 
As N — > oo, and g > 1 zs /ized, 



2g+l 



(5.127) PI (i>*,i^) e ( c , d)J ~ (2TT)- — ^/T q 
where the matrix @*(q) is defined as above, and 



exp 



(c,d) 



z T e*(q)^ )d~z, 



(5.128) 



T q = T{p + 2) m 



Proof: We follow the known technique for passing from a local limit theorem to an integral one, that 
is exposed in detail in [37], p. 60 (see also |10j . p. 80 and [12]). Summing (|5.124j) over (x°, y°) 6 Gat, 
we obtain 



y*,K*) 6 (c, d 



(27T)" 



2g+l 



/o(p-l) 



So 



1(9+1) 



/o(p+l) 

^0 



(5.129) 



cxp( --^ T G*(g)^) (l + e(^,JV)) 



iV -> oo. 



From the preceding asymptotic formulae we derive the crucial fact that in (|5.129p . 



(5.130) 



sup e(~z, N) 0, iV oo. 

"5*e("c , d ) 



Based on this property of uniform convergence, we treat the RHS of (|5.129p as a Riemann sum with 



the asymptotically equidistant spacings I \G 



,0) 

r N 



/o(p-l) 
So 



5%, j = 0, . . . , q and ( \G 



'N 



. p+2 



i = g + l,...,2 g! asiV^oo. 
Next, we have to show that the main term in the RHS of (|5.129p is indeed the gaussian distribution. 
This is equivalent to proving that 



(5.131) 



dete*(g) =T g , q>l, 
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where T q is as in f|5. 128j) . By the definition (|5. 106|) of Wj, 



(5.132) 



+ 



ff(p) 



/i(P + l) 



/i(P"l) ' //(P-IM /;(P-1K 

We first prove (|5.13ip for g = 1. Using (|5.132[) we have from (|5.124p 



j = 0, . . . , q. 



(5.133) 



e*(i) 



/ /o(p+1) 
/o(p-l)w)() 



MP) 
\ /o(p-l)w) 





/i(p-l)u>i 

/iCp) 

/i(p-l)wii 



/o(p) 
/o(p-l)w)o 

flip) 
/i(p-l)tui 

U)Q 11)1 



Now the claim for g = 1 follows from the identity /o(p+l) + /i(p + l) = CT(p + 2), p > and some 
algebra. This shows that the random vector (z>q, £*, X*) weakly converges to the gaussian random 
vector with zero mean and the covariance matrix (0*(1)) _1 . This implies (see [TU], p. 89, Theorem 
2.6) that for q = 1 the sequence of the distribution functions of the random vector (j>q , i>\, K*) is tight 
under each u\ > 0. Consequently, by the definition of tightness and the form of the stratifications 
considered, we deduce the tightness of the distribution functions of the random vector (z>*, K*), for 
all q > 1. Finally, by Prohorov's theorem, the limiting distribution of the above vector is the one of 
a probability measure. This proves (|5.13ip for q > 1. 

Now we derive the explicit form of the covariance matrix in Theorem 4.1. For a given q > 1, let 
(Q*(q)) 1 = {^mfc(^)}m fc=o ^ e ^ e covariance matrix of the limiting gaussian distribution in (|5.127p . 



Lemma 5.7 



r( P +2) • 



(5.134) r mk (g) = { 



fm(,P l)lm= k 



f ( n \1 , _ f m {p)fk- q (p+l) 
Jra\P) s -m=k-q r(p+2) : 



fk{p)lk=m-q 
Jm-q(P + l)lm=fc 



fk(p)f m - q (p+l) 

r( P +2) 

/m-g(p+l)/fc-,(p+l) 



i/ < m, k < q, 
if < to < q, q - 
if < k < q, q 4 



- 1 < fc < 2q, 
1 < to < 2q, 



r( P +2) 



i/ q + 1 < m, fc < 2g. 



Proof: To understand the technique behind the inversion of the matrix 0*(g), it is convenient to 
verify (|5.134p first for q = 2. This is easy to do with the help of the identity q - =G fj (p + 1) = 
CT(p + 2), which holds for any q > 1. Taking into account the structure of the matrix 0*(g), the 
verification for q > 2 can be done in the same way. ■ 

This completes the proof of the statement (ii) of Theorem 4.1.. We note that by (I5.134p . all covariances 
^mfc(<z) between the scaled component counts Uj, as well as between masses are negative. 

The convergence of moments, as stated in (iii) of Theorem 4.1 results from (|5.130p . For the proof, 
it is needed only to replace probabilities with moments in the argument given in [37], p. 61 (see also 
[23], p. 67). We see from (|5. 134|) that for the proof of the claim (iv) of Theorem 4.1 one has to 
examine the behavior in k — > oo of the integrals fk(p), fk-q(p+ 1) in the case when Uj, j = 0, . . . ,q 
are equidistant points and q — > oo. Now Theorem 4.1 is completely proved. 
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We proceed to the proof of Theorem 4.3. First, we recall the scaled quantities Kj in (|4.45p 
and denote 



(5.135) 



K = (K 1 , 



It follows from the definition of ua and Kj, j = 0, . . . ,q that 

<? q 



(5.136) 



y 0< 



V 3 



By Theorem 4.1, ()5.136p and (|5. 134|) . the random vector ( £> , K) converges weakly to the gaussian 
vector with zero mean and the covariance matrix T(q) = {e m k}^ k=0 given by: 



(5.137) 



6 rn k 



' fjy 11 f " (p ^) 

J0\P *-) T(p+2) ■ 
Mp)b P {u k ) 

i>+2) ' 

fo{p)bp+l{u k - q ) 

r( P +2) 

^p— 1 (^max(m,/c) ) 
^p(^max(mjfe- q) ) 



k ^p+1 ("maxfmi)-?) 



bp{u m )b p {u k ) 

T(p+2) ' 
b p (u m )b p+1 (u k _ q ) 
r(p+2) ' 

_ fcp+lOm-g)frp+lOfc-<,) 

r( P +2) ' 



if m = k = 0, 

if m — 0, < fc < <j, 

if m = 0, g + 1 < k < 2q, 

if < m, k < 

if < m < q, q + 1 < k < 2q, 

if q + 1 < m, fc < 2g. 



Applying now (|5.13Tj) with < m,k < q gives the claim (i) of Theorem 4.3. 

For the proof of the claim (ii) of the theorem, we first recall that in the notation (|2.1|) . Uj = 
v(rNUj), j = 1, . . . ,q, with rjy as in (|4.43|) . Thus, letting q = 1 and u\ = u, we rewrite the first 
part of (|4.45p as 

v{r N u) - St 



(5.138) 



V/o(p-l) 



to get with the help of ([SMD, ([EMD, (|Q31 and the notation the desired relation (^35]) between 
zv'j and the fluctuations A(«j) for j = l,...,q. By virtue of the convergence stated in (i) of Theorem 
4.3, the fact that NrJ^ 1 — > oo, as iV — > cxd and the continuity of Varu\ in u, we derive that v\ is 
stochastically equicontinuous (see e.g. [33]) p. 471) in u € [c, d\, where [c,d] is any finite interval. 
This implies the claim (iii) of the theorem. 

Now we turn to the proof of Theorem 4.6 which will be done in the manner similar to the proof 
of Theorem 4.1. For a given q > 1, we consider the stratification points Mj,j = 0, . . . , q + 1 as in 
(|4.56p . and define the quantities S* , Ej ,V* , j = 0, . . . , q, as in (|5.66p and (|4.47j) . Clearly, Lemma 
5.1 remains valid with the above points Mj. In accordance with the problem considered, the vectors 
L* and N* are taken now in a slightly different form: 



(5.139) 
and 

(5.140) 



Vj = Sj + X2jJSj, j 



0, 



N 
j 



Ej + x 2 j+iJV 



0,...,q-l, N* = N 



g-l 
3=0 



Limit shapes 



33 



where xq, ■ ■ ■ ,%2q are arbitrary reals not depending on N. Under 5 = 5n chosen as in (|5.9ip . the 
following analog of Lemma 5.3 is then valid: 



(5.141) 



R(L*,N* 



W 



,+i mi / — i ex p i 



x 



2j ' 



N oo, 



(2vr) - \,. M x /.^ / \, j=0 
where is defined by (|5.79p . It is important to note that here, in contrast to the previous setting, 



(5.142) 



lim SMj=0, j = 0,...,q, 



by (|5.9ip and (|4.56p . Next we will intend to show that the local limit theorems (|5,85p are valid with 
Kj replaced with r? given by (15.1490 below. Firstly, by the definition (14.561) of Mj and by (I5.142p . 
the integral test gives for r > and N — > oo, 



(5.143) Yj a k k 

k=Mj 

From this and (|5.65p we derive that, under N —* oo and r > 



r(p + r)5-P" r , if j 



if j = 0,...,q-l 



(5.144) 



0(Mj +1 ), if j = 0,...,(?-l, 
0(<r r ), if i = g- 

In the case considered the analog of (|5.97p will be 



to,j - ' 



(5.145) 

Thus, we obtain that for any r > and N — ► oo, 
(5.146) (E/5J) t 0J 



L*J 2 Mr + \ log 2 £* if j = 0, . . . , g - 1 
6(L*yhog 2 L*, tij = q. 



i O[ M j+ l log 2r M j+1 ) , if j = 0, . . . , q - 1, 



O ( <5lr log 2r <5 



if j = q. 



It is easy to see that (I5.146P guarantees the conditions (|5.10Q[) and (|5. 101[> . Next, we have from 
(153431 



E* 

(5.147) lim - 3 == 



lim^eo (*@Pp_) ^ ±i= ^ := if j = 0,...,g-l, 



r( P +i) 



Vr(p)r( P +2) 



if J = 3, 



where the limits < a, < 1, j = 0, . . . , q, by (|4.56p and some algebra. It follows from (|5,140p that 



(5.148) 



j=0 
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We obtain from (loTTO . (I436D . (loTTiOl) and (I5TT391) 



iv; - EyCi) _ at; - -girL* 



I TP* 

x 2 j+i, V* - x 2j 



X2j+1 



x 2j 



X2j+1 ~ OtjX 2 j 



1-af 



(5.149) 



jV^oo, i = 0, 



while for j = q we have 



v / Vary(9) 



g-l 

i=o 

v / Vary(<?) 



Z2g 



(5.150) 



S**Vary(9) 



g q x 2q 

1 - 



S£Vary(<?) 



T q , N — > OO. 



We observe that \Tj\ < oo, j = 0, . . . ,q. Proceeding further along the same lines as in the proof of 
Lemma 5.4 leads to the desired analog of (|5.85p . Consequently, (|5.122j) takes the form 



(5.151) fl (Lvyr*) ~ | n 



exp 



E 

j=0 



(2tt)^ \ j=0 y/StVarYV 
Now we define the scaled quantities for the problem considered 



2,2 
X 2 j + Tj 



N — > oo. 



K*-E* 



j = 0, . . • , q, 
j=0,...,q-l 



and denote 



*>* = (/>*,...,*>*), K* = (K*,...,K*_ 1 ). 



(|5.147D and (15.121j) imply that 



VarZ N ~ V* and 



(5.152) 



Va^^l^b <0 °' JV - >00 ' J = °'--- 
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Substituting (j5.152j) into (|5.15ip gives 

" = ~x, k 



y 



2(7 + 1 

(2vt)- — 



(5.153) 



exp 



£ 

3=0 



2 . (^2j+l - OtjX 2 j) 



3=0 
2\ 



n 



spy; \f=oJ^ 



x 2j + 



I -a] 



+ 



2g 



1 - «2 



iV -» oo, 



where we set 1? = (xq, »2, • • • , ^2?) and y = (xi, x$, . . . , x 2q -±). 



Remark 5.8 At the first glance, the expression \5.153\) and its analog \5.t2J$ for Theorem 4-1 look 
very much alike. A crucial difference between them is that in \5.12J$ the variables x° are linked via 

2q+l 

the relation x° = 0, while in ( 5.124) the variables Xj are free. Formally, this is implied by the 



i=9+i 

fact that the value K q in \5.101^ depends on x° and x-L+n whereas the value r q in 15.150]) depends 
on x 2q only. As a result of the aforementioned difference, conditional independence in Theorem 4-1 
transforms to independence in Theorem 4-3. 

From (|5.153p . the claim (i) of Theorem 4.6 follows immediately. It is left to find explicitly the 
covariance matrix of the corresponding limiting gaussian distribution. For the sake of convenience, 
we write z = (xq, . . . , x 2q ). Then the matrix, say B, of the quadratic form in the exponent of (|5.153p 
has a diagonal block structure with the blocks Bj of the form 

l «. 



l-ai l-aj 



and 



(5.154) B q - — 



\ — a 2 ' 



Thus, 

(5.155) det B = f[ det Bj = f[ - 1 ^ 



2 

3=0 3=0 * 3 



and 



(5.156) Bj 1 = ( 1 a > Y j = 0, . . . , q - 1, B~ x = 1 - a\. 

In accordance with the above notation for z, we set x = (Aj> ^o> • • • > ^q— i> ^q-ii ^q) w hich 

enables to rewrite (|5.153p as 



p ( * = f ) — ^= n -7= -p (4-* b ^ t ) ■ 

(5.157) iV ^ oo. 



3 3 
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Hence, the covariance matrix B 1 of the gaussian distribution can be written as the Kronecker 
product of blocks BJ as in (|5.156p : 

(5.158) B~ l = Bq 1 ® Bl 1 ® . . . ® B~\ 

Finally, by the same argument as before, we derive the weak convergence as stated in (ii) of Theorem 
4.6. 

6 Comparison with limit shapes for the generalized Bose-Einstein 
and Fermi-Dirac models of ideal gas 

Recall (see Section 3) that the generalized BE and the FD models in the title are associated with 
the generating functions T^ BE ^ , T^ FD ^ respectively: 

fc>l V ' k>l 

(6.159) \x\<l, m k = Ck p ~ 1 , p>l, C > 0. 

Limit shapes for these models, say C^_f ,C£^ respectively, were obtained by Vershik in 



p—cix roc p-cix 

(6.160) C£?\u)=l x^——dx, C^\u)= x^——dx, u>0, P >l, 



where c\ = c\{p) > 0, C2 = ozip) > are normalizing constants. 

It is interesting to compare these limit shapes with the ones, denoted l p -i, in our Theorem 4.3: 

f 00 e~ x x^~^ dx 

(6.161) lp - l{u)= T(p+1) ' P> °- 

Firstly, note that the limit shapes for both models are derived under the scalings of the same 
i 

order NP+ 1 . In view of Proposition 2.5, the immediate explanation of this striking coincidence is that 

by the result of Vershik and Yakubovich (|47j, Section 5), the models (|6.159p of ideal gas have the 
i 

same threshold NP+ 1 as the models considered in the paper. In a broader context, we point that the 
generalized BE models and the models (|2.9p studied in our paper, are linked with each other via the 
exponentiation of the generating function !f^ BE ^ (see [5]). Moreover, the Bell-Burris Lemma 5.1 in 
[7] states that if the parameter function m = {m^} for the generalized BE model (multiset) is such 
that lim/j^oo J^ 1 = y with y < 1, then the corresponding multiplicative measure for the BE model 
is asymptotically equivalent to the measure (assembly) given by (|2.9p with ~ m^, k — > oo. 
However, for the generalized BE models (|6.159p . lim/^oo m £ i ~ 1 = 1, which explains the following 
basic difference in the form of the limit curves C^_f ; and Zp_i. From (IBTTBTT) . (I6TT60|) we see that 

(6.162) /p-i(0) =p~ x < oo, p>0, 
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while 



(6.163) C^ BE) (0)= lim -ll]og(l-exp(-^X) = 00, C p _i(0) < 00, p > 1. 

u^0+ IT \ V v 6 / / 



By (|2.ip . (|2.2p and (|2.3p . the value of a limit shape at u = "approximates" the random variable 
z?(0) = the total number of components in a random partition, multiplied by the factor In [44] . |43| 
the phenomenon (|6.163p is linked with the Bose -Einstein condensation of energy. According to this 
interpretation, the finiteness of the limit shape at u = indicates the appearance of condensation of 
energy (around the value ^-), and in view of this the value p = 1 (=uniform distribution on the set 
fijv) was distinguished as the phase transition point for the generalized BE models ()6.159p . We now 
give an analytic explanation of the fact that Cq (0) = 00 in contrast to £q(0) < 00. By the seminal 
result of Erdos and Lehner(1941) (for references see [47J), the number of components in a random 
partition of N is asymptotically ^y/Nlog N S> y/N = On the other hand, it was proven in 
that the number of components in our model with p = 1 is asymptotically y/~N = h-^-, where 
constant h is as in (|4.43|) . 
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